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Abstract — We Investigate a wireless network localization sce- 
nario in which the need for synchronized nodes is avoided. It 
consists of a set of fixed anchor nodes transmitting according 
to a given sequence and a self-localizing receiver node. The 
setup can accommodate additional nodes with unknown positions 
participating in the sequence. We propose a localization method 
which is robust with respect to uncertainty of the anchor 
positions and other system parameters. Further, we investigate 
the Cramer-Rao bound for the considered problem and show 
through numerical simulations that the proposed method attains 
the bound. 

Index Terms — Wireless sensor networks, maximum a posteriori 
estimators, Cramer-Rao bound, anchor uncertainty. 



I. Introduction 

The rapid emergence of wireless sensor network (WSN) 
applications, often requiring accurate knowledge of the node 
positions, has created remarkable research interest for the 
WSN localization field |1|. In such applications, nodes need 
to estimate their own position, i.e., perform self-localization, 
by processing noisy range measurements with respect to a set 
of anchors in a decentralized fashion. 

Typically, the anchor nodes positions are assumed to be 
known exactly, thus neglecting a potential source of error 
in many application scenarios. In the literature, the problem 
of sensor network localization in the presence of anchor 
position uncertainty has been posed as a maximum likelihood 
(ML) estimation problem, cf. [2|. However, due to the highly 
nonlinear likelihood function, a closed-form expression of 
the ML estimator cannot be derived. Therefore the solution 
has been approached using semi-definite programming in [2 1 
and second-order cone programming in Q. Further, in (3], 
the authors studied the joint network localization and time 
synchronization problem with inaccurate anchors, pointing out 
the strong interdependence between the timing and positioning 
aspects. All the above mentioned works assume that the timing 
measurement noise variance is known. 

In this letter, we consider the problem of self-localization in 
sensor networks with uncertainty in anchor position, without 
assuming knowledge of the timing noise level. More impor- 
tantly, we propose a system configuration which allows for a 
receiver node to perform self-localization even without time 
synchronization of the anchor nodes. In such a configuration 
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the transceiving anchor nodes transmit signals in a prede- 
termined sequence, one after the other, with a turn-around 
delay which we assume is not perfectly known at the receiver 
By passively listening to the transmissions and exploiting 
available prior knowledge, the receiver can localize itself as 
well as the transceivers that participate in the sequence. This 
setup can accommodate the self-localization of an indefinite 
number of receiver nodes. 

We propose an iterative maximum a posteriori (MAP) esti- 
mator that effectively copes with uncertainty in the deployment 
of anchors and with tolerance of the turn-around delay of 
the asynchronous transceivers. It also provides delay estimates 
which may be useful for hardware calibration. Moreover, we 
study the fundamental performance bounds for the problem 
considered by deriving the hybrid Cramer-Rao bound. Finally, 
we evaluate the performance of the proposed estimator by 
numerical simulations of an ultra-wideband sensor network 
setup as an application example |5|. 

Notation: ||x||-w = Vx^Wx is the weighted norm, where 
W is positive definite. 

II. Problem formulation 

We consider a wireless network of A — 1 transceiving nodes 
operating asynchronously with local clocks. They transmit 
signals according to a known sequence, denoted T, set across 
the network |6|; when the next node in the sequence receives 
a signal, it transmits in return after a certain delay. On this 
basis, the goal is to achieve self-localization of the Ath node, 
which is a passive receiver that knows T. 

The signals are assumed to have a resolvable temporal 
signature that allows for timing events, e.g., pulses, symbol 
boundaries, etc. and the propagation velocity c is known. Let 
Xj € M'' denote the position of node i, where d = 2 or 3, 



and 



|xi — Xj||2 denote the range between nodes i and 



j. Then the observed time interval between a pair of signals 
received at node N, involving transceiving nodes i and j, is 
modeled by 

1 „ 1 1 



y(*'j) = -Pij + Sj + -pj,N Pi.N 

c c c 
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where Sj denotes the turn-around delay at node j, cf. (6\, fT\. 
The delay dj is generated without any common time reference 
between the nodes, and is therefore asynchronous. The noise 
y;(' j) arises from three uncorrected timing measurements, one 
at node j and two at node N, and is modeled as zero-mean 



Gaussian with unknown variance E[(w(*'^^'i^ 
observed time interval 



] = a . The next 
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uses one timing measurement from the previous observation. 
Modeling the timing noise variance across nodes equally we 
have the correlation E[i(;(''^^u'(-'^'^^] = cr^/3 for all consecutive 
observations. 

Prior knowledge of the positions will be modeled as x^ ~ 
A/^(/i.j, Pi), where /x, and Pi are known, \/i. Setting P^ — 0, 
leads to a noninformative prior, p(xi) ex 1 |]8l. Further, to 
avoid signal collisions it is necessary that the delays exceed 
Pmax/c, where pmax IS the maximum range between any pair 
of transceivers and can easily be ensured in any bounded 
localization scenario. Whilst the delay in each transceiver node 
may be set to some nominal value /i^-, the actual delay 5j 
wiU deviate due to hardware imperfections. We model this 
as 5j ~ NIhstCtI), with a given a^. For the unknown noise 
variance, we assume a noninformative prior p(cr^) oc 1/ct^ |8|. 



For notational simplicity we write 9 = [xj^ • • • x 



idN 



N\ 



and 6 =[5i- ■ ■ 5n-i] ' G R^ ^- The goal is to estimate 6, S 
and a^ from a set of M observations {ym}m=i- 

III. MAP ESTIMATOR 

For a given sequence T, the pair of nodes involved in each 
observed time interval {ym}m=i is known. In vector form, the 
observation model is 

y==c-iHg(i?)+weM^^ (3) 

where i? = [6^ J^]^ e M^, and T = diV + iV- 1. The noise 
follows w -- 7V(0,CT^Q), where [Q]iA = 1, [Q]i,j = ^, Vi, j 
such that |i — j| = 1 and [Q],;.j = otherwise. The nonlinear 
mapping g(i5) = [p'^ (9) <5^]T e ]^n{n-i)/2+n-i contains 
all the pairwise ranges and delays, and H is determined by 
the transmission sequence T, cf. ([U. We aim to find the 
maximum a posteriori (MAP) estimator, i.e., the maximizer 

A. Concentrated cost function 

Using Bayes' Rule, the MAP estimate can be computed by 
maximizing J{'d,a^) — \np{y\'d,a^) + Inp(cr^) +\np{'d). 
Further, define 



M + 2 



1 



B. Iterative solution 

The optimization problem (|6]l does not lend itself to a 
closed-form solution. A standard method for solving (|6) is gra- 
dient descent. However, its performance is heavily dependent 
on the user-defined step length. Instead we propose an iterative 
solution by first exploiting the linearization around an initial 
estimate '&i, i.e., g(i9) ~ g{'di)+T{'di)&, where ^ = -d-hi 
is the iteration increment and T{'d) is the Jacobian of g(i9). 
Then we may write a cost function which approximates (|7]i as 



where Ki is a constant. Similarly, 



2 

p-1 



K^ 



(5) 



where K2 is a constant. Maximizing ^ with respect to cr^ 
yields the estimate ct^ = ||y - c"iHg('j9)||2 _j/(M + 2). 
Inserting this back into (|4]i, and combining with Q, results 
in a concentrated cost function. The MAP estimator is then 
given by 

i^ = argminy(i9), (6) 



^\\^,-'^\\l-^, (7) 



where 

y(,9)^iln||y-c-iHg(,9)|pQ_ 

/3 = 1/(M + 2), ^l=[^ll ■■■ fij, fislJj_,V and P^ 
diag(Pr\---,P^\a7'ljv-i). 
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where we introduce g^ = g('i?£), T^ = r(t?£), yi = 
y — c^^Hgc, Ge = c^^HTi and fig = fi — t?^. We aim 
at iteratively updating the estimate -de by finding the optimal 
increment i9. 

To this end, we use a fixed-point iteration. First, we note 
that the gradient of V^ (1?) equals 

d^Vi = a{^)GjCl-\G,d - ye) + (3P-\-d ~ fi,), 

where a(^) = l/||y<> - Gf^||^_i. Then, we hold a{d) fixed 
and solve di}Vi — Q, thus obtaining the following fixed-point 
iteration 



i9 



= {a[^)Gjq,-^Gt + l3P- 
x(aWG/Q-iy, + /3P-iA, 



(9) 



By iteratively applying (|9]l, starting with the zero increment 
1? = 0, we converge to a stationary point. The analytical 
convergence properties are difficult to derive and are beyond 
the scope of this letter However, in Section |V] we show the 
convergence properties by numerical evaluation in a practical 
scenario. Further, we provide a practical method for obtaining 
an initial estimate 1^0 • The iterative estimator is summarized 
in Algorithm [T] where e is the convergence threshold. 

Algorithm 1 Iterative MAP estimator 



2 ln^'-^lly-c-^Hg(i5)||2j_,+Xi, 

(4) 



1 


Input: y, /x, P, c and i9o 


2 


Set £:=0 and;3= l/(A// + 2) 


3 


repeat 


4 


y^ = y - c-iHg(i9^) 


5 


Gf = c-^nvi^t) 


6 


fif = fj.-'de 


7 


i9 :=0 


8 


Repeat ^ until convergence 


9 


,9^+1 =-d + 'de,£:^£ + l 


10 


until W-de — i?£-i 2 < £ 


11 


Output: ■&£ 



IV. Cramer-Rao bound 

e M-^+^ and r) be any estimator 



Let T7 ^ [0^ S^ cr^ 



that is conditionally unbiased with respect to the deterministic 
parameters. Then its mean square error (MSB) matrix is 
constrained by the hybrid Cramer-Rao bound (HCRB) |^, 
Cij >z J^\ where J^ = J^ + J^ e 



d(T+1)x(T+1) 



Here J^ = E^[J^(?7)] is the expected Fisher information 
matrix, where f) denotes the subset of parameters that are 
modeled as random quantities and 



[J^(r/)].,, = 



H'Q^H— ^ 



(T^ d-qi drjj 2cr* dT]i drjj 

as given in |10J. As the expectation does not have a closed 
form solution, we evaluate it by Monte Carlo simulation. If 
a subset of node positions, 0„, and the noise level, a^, are 
treated as deterministic and unknown parameters, while the 
remaining parameters, 9a and 6, are random Gaussian then 
the prior information matrix is given by 



0„ 























cr^^ljv-i 












n(T+l)x(T+l) 



where P^^ is the covariance matrix of 9a- This division 
between deterministic and random parameters may be useful to 
study practical configurations where we lack prior knowledge 
on the position of a subset of nodes. 

V. Numerical results 

To evaluate the performance of the MAP estimator proposed 
in Section |lll] we provide numerical simulation results in 
an ultra wideband wireless sensor network scenario [5]. In 
particular, we compare the HCRB derived in Section ITVl with 
the root mean square error (RMSE) of the MAP estimates. 

A. Setup 

We consider a network of A^ = 6 nodes consisting of 
one self-localizing passive receiver, without prior knowledge 
of its position, and Na = 4 anchors with Pg^ = aaT-2Na- 
In addition, we consider one auxiliary node, which is a 
transceiver participating in the sequence with a noninformative 
position prior Therefore, the number of unknown-position 
nodes is Nu = 2. The simulation has been performed using 
the topology shown in Fig. [T] and setting fis = 10^^ s. 
In each realization, the positions of the anchors and the 
delays have been randomly generated according to their prior 
distributions. Furthermore, we construct the sequence T to 
ensure that all pairwise combinations appear at least once. 
E.g.,r={l, 2, 1,3, 1,4, 1,5,2,3...}. 

The average RMSE of the position and delay estimates are 
given by 



RMSE^ 



N^ 



y^{-^Y 



where ^ can be either 9u, in which case N^^ = Nu, or d, and 
N^ ^ N -I. Here C| is the MSE mati'ix of ^. We estimated 
the RMSE from 10^ Monte Carlo iterations. 

We initialize the MAP estimator of Algorithm [T] with 



i9n 
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where u is 



an arbitrary position initialization of the auxiliary node, with 
u 7^ /Xj. The self-localizing node's initial position estimate p, 
is the centroid of the anchor positions, i.e. fi — ^ 12i=i f^i- 



Further, we set e = 10 
Step 8 in Algorithm [T] 



The same tolerance is used for 
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Fig. 1. Scatter plots of the true node positions and error ellipses for all 
nodes estimated at the self-localizing node indicated by (*). The anchors are 
denoted by (x), the auxiliary node by (A) and the anchors centroid by (D). 
The solid black ellipses indicate the CRB and the dashed red ellipses indicate 
the MSE performance of the MAP estimator For visual clarity the sizes have 
been scaled to correspond to 99% confidence ellipses of a zero-mean Gaussian 
distribution. Here a = 2 ns, (Ja = 0.2 m and ag = 10 ns. 
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Fig. 2. Average RMSE of the position estimates of all unknown-position 
nodes at the self-localizing node as a function of the standard deviation of 
the measurement noise a. Two values of the prior on the anchors positions 
are shown. Here, erg = 10 ns. 



B. Results 

Fig. [T] shows the topology of the network, together with the 
error ellipses produced by the MAP estimator and given by 
the HCRB. It can be seen that the proposed estimator attains 
the bound. Next, we provide an evaluation of the position 
estimates as a function of timing noise level a. In Fig. |2] 
we vary the uncertainty of the anchor positions which is 
parameterized by Ua. In a low noise scenario, we achieve an 
accuracy of the same order of magnitude as the anchor position 
prior. In Fig. |3] we vary the uncertainty of the delays set 
by as- The results show that the proposed position-estimation 
method is robust with respect to deviations from the nominal 
delay. Similarly, Fig. |4] shows the robustness of the delay 
estimator We tested a range of deviations that spans two 
orders of magnitude. Furthermore, the extended simulation 
setup in Fig.|5]shows that the proposed method can localize 10 
auxiliary nodes which participate in the transmission sequence 
and are not collocated. However, the MAP does not attain the 
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Fig. 3. Average RMSE of the position estimates of all unknown-position 
nodes at the self-localizing node as a function of the standard deviation of the 
measurement noise a. Two values of the prior of the delay are shown. Here, 
(Ta = 3 cm. 



Fig. 5. Scatter plots of the true node positions and en'or ellipses for a uniform 
random placement of 10 auxiliary nodes denoted by (A). Here cr = 2 ns, 
Oa = 0.2 m and a^ = 10 ns. 
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Fig. 4. Average RMSE of the delay as a function of the standard deviation 
of the measurement noise a. Two values of the prior of the delay are shown. 
Here, Ua = 3 cm. 



HCRB in all scenarios, in particular when nodes are close 
to anchors as can be seen in Fig. |5] Moreover, Fig. |6] shows 
the convergence behavior of the estimator, the average is 5.12 
iterations. 

Finally, we repeated the simulations in the scenario in Fig.lT] 
assigning a more conservative prior with respect to the actual 
variability of the anchor positions. In particular, we defined 
(7^ ~ lOcTa, where a'^ is used by Algorithm [T] and Oq, — 0.2 m 
is used to generate the random anchor positions. We obtained 
RMSEe„ = 0.37 m and HCRBg^ = 0.26 m, indicating 
inherent robustness to such model mismatches. 

VI. Conclusion 

In this letter, we have considered a wireless network lo- 
calization scenario with asynchronous nodes where the entire 
network topology can be estimated by a self-localizing receiver 
node. We have proposed a MAP estimator and compared its 
performance with the HCRB that we derived for the problem. 
The simulation results show that the estimator attains the 
bound and is robust with respect to uncertainty of the anchor 
positions, delays and noise level. In addition to the anchor 
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Fig. 6. Histogram of the number of iterations of the MAP estimator in the 
scenario in Fig. [T] for 10'^ realizations. 



nodes, we found that the setup allows for the localization of 
auxiliary nodes participating in the transmission sequence. 
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